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ABSTRACT 

Vibrations generated by machine driveline components can cause excessive noise and structural dam- 
age. Magnetostrictive materials, including Galfenol (iron-gallium alloys) and Terfenol-D (terbium-iron- 
dysprosium alloys), are able to convert mechanical energy to magnetic energy. A magnetostrictive 
vibration ring is proposed, which generates electrical energy and dampens vibration, when installed in 
a machine driveline. A 2D axisymmetric finite element (FE) model incorporating magnetic, mechanical, 
and electrical dynamics is constructed in COMSOL Multiphysics. Based on the model, a parametric 
study considering magnetostrictive material geometry, pickup coil size, bias magnet strength, flux path 
design, and electrical load is conducted to maximize loss factor and average electrical output power. 
By connecting various resistive loads to the pickup coil, the maximum loss factors for Galfenol and 
Terfenol-D due to electrical energy loss are identified as 0.14 and 0.34, respectively. The maximum av- 
erage electrical output power for Galfenol and Terfenol-D is 0.21 W and 0.58 W, respectively. The loss 
factors for Galfenol and Terfenol-D are increased to 0.59 and 1.83, respectively, by using an L-C resonant 
circuit. 
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1. INTRODUCTION 

Vibrations generated by machine driveline components can cause excessive noise and structural damage. 
Piezoelectric materials are able to transduce mechanical energy to electrical energy; thus they can be 
implemented as passive structural dampers. A PZT device was proposed which applies damping in a 
driveline to reduce unwanted vibration and noise. 1 The device has the appearance of a metal ring, as 
shown in Figure 1. 

Magnetostrictive materials, which couple mechanical and magnetic energy, have been widely imple- 
mented in damping and energy harvesting applications. 2 ~ 7 Many magnetic materials exhibit magne- 
tomechanical coupling, but only a few materials containing rare earth elements show significant or giant 
magnetomechanical coupling effects. Terbium-iron-dysprosium alloys, known as Terfenol-D, are able to 
generate 1600 ppm magnetostriction, but a high magnetic field (160 kA/m) is required to saturate this 
material. 8 Iron-gallium alloys, known as Galfenol, are a recent class of magnetostrictive materials that 
exhibit moderate magnetostriction (400 ppm) and high magnetization (1200 kA/m) at low magnetic 
fields. 9 Figure 2 shows a new vibration ring based on magnetostrictive rods. The rods are arranged 
radially around the circumference, and are made of Terfenol-D or Galfenol. With this configuration, the 
vibration ring is able to convert radial vibration energy to electrical energy, thus creating a vibration 
damping effect. Unlike Terfenol-D and piezoelectric materials, Galfenol has a mechanical robustness 
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Magnet 


Figure 2. Compression spoke design of a magnetostrictive vibration ring. Terfenol-D has a tensile strenght of 
about 28 MPa, thus it must be operated in compression. The Terfenol-D rod is press fit between the two rings. 
Galfenol can be threaded into the rings, since it has a tensile strength of about 500 MPa and good resistance to 
torsion and impact. 


Proc. of SPIE Vol. 9433 94330C-2 


Downloaded From: http://proceedings.spiedigitallibrary.org/ on 04/07/2015 Terms of Use: http://spiedl.org/terms 




Table 1. Comparison of Terfenol-D and Galfenol. 6 ’ 12 14 


Material 

Terfenol-D 

Galfenol 

Composition 

Tb x Dyi_ x Fe 2 (x « 0.3) 

Fei_ x Ga x (x = 0.12 - 0.33) 

Conductivity (S/m) 

1.72 x 10 6 

5.96 x 10 6 

Relative permeability 

3 - 10 

1 - 700 

Saturated Young’s modulus (GPa) 

115 

60 - 80 

Density (kg/m 3 ) 

9250 

7870 

Tensile strength (MPa) 

28 

500 

Static hysteresis loss per cycle (kJ/m 3 ) 

23 

2.1 (Constant field biasing) 


similar to steel, and thus can be machined, welded, and formed. Without undergoing permanent depo- 
larization, magnetostrictive materials also maintain significant magnetomechanical coupling over a broad 
temperature range. 10,11 Table 1 compares the mechanical and magnetic properties of Terfenol-D and 
Galfenol. 
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Figure 3. (a) Geometry in COMSOL Multiphysics for Galfenol analysis; (b) geometry in COMSOL Multiphysics 
for Terfenol-D analysis. 


Instead of modeling the entire vibration ring directly, this study investigates the damping and energy 
harvesting properties of the Galfenol and Terfenol-D rods, which are the building blocks of the design. 
Two dimensional (2D) axisymmetric models, as shown in Figure 3, were created in COMSOL Multi- 
physics. The length of the rods is 10 mm. The top and bottom platens are made of stainless steel A430. 
The return path is made of 1018 low carbon steel. A 750 Hz, 280 N amplitude sinusoidal force (the 
target vibration source) is applied on the top surface of the rods, while their bottom surface is fixed. 
Cylindrical Alnico magnets (Grade 8) with a remanent flux density of 0.6 T are used to generate a bias 
magnetic field through the magnetostrictive rods. The 2D axisymmetric model incorporates magnetic, 
mechanical, and electrical dynamics. Based on the proposed model, the performance of the Galfenol 
and Terfenol-D rods is investigated in terms of loss factor due to electrical energy loss rj and average 
electrical output power P. A parametric study for rod geometries, coil size, and magnetic flux path is 
also presented. 


2. PERFORMANCE CONSIDERATIONS 


2.1 Loss factor r/ 


Loss factor is a common measure of a material’s intrinsic damping, 
oscillation as 15 

D 

V ~ 2t tW’ 


It is defined for a steady-state 


( 1 ) 
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where D is the energy dissipated and W is the energy that is returned to the input with each cycle. For 
magnetostrictive components, there are several loss mechanisms, including hysteresis loss, mechanical 
damping, eddy current loss, and dissipation of generated electrical energy. In this study, the only loss 
considered is the latter. Therefore, the electrical energy loss per cycle, £) oss , is maximized in order to 
maximize the loss factor. 

In the frequency domain, the sinusoidal input force F* with a frequency /o can be written as 

F* = F 0 e^ F , (2) 


where Fq is the force amplitude and <pp is the phase of the signal. Since magnetostrictive materials are 
highly nonlinear, the tip displacement due to F* consists of a fundamental component S* and higher 
harmonics. The loss factor in this study is calculated for the fundamental frequency of the response 
using the Fast Fourier Transform (FFT) method. 16 The fundamental component of the displacement is 
defined as 

S* = S 0 e j4,s , (3) 

where Sq and <ps are the displacement amplitude and phase, respectively. Following ASTM D5992- 
96(2011), 16 the loss factor is calculated from the phase difference as 

V = tan(^ F - <j> s ). (4) 


The ASTM standard method for calculation loss factor is employed so that the results are comparable 
with commonly available data on passive material damping, such as for rubber or metals. 


2.2 Average electrical output power P 

Figure 4 shows the equivalent circuit for the electrical system, including the resistance and inductance 
of the pickup coil and electrical shunt components. Both purely resistive and parallel resistive-capacitive 
shunts are analyzed in this study. The average electrical output power P on the resistor R ^ quantifies 
the harvested energy and it is calculated as 


P = 



Vi M‘ 

Rl 


dt , 


(5) 


where Vi(t) is the voltage on the resistive load and T is the period of the vibration. 
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Figure 4. (a) Resistive load; (b) capacitive load. 
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2.3 Purely resistive load 

A previous study 1 ' showed that P for Figure 4(a) is maximum when 

R l = y/B? c + (M) 2 , (6) 


where R c is the coil resistance, L c is the coil inductance, and w is the radial frequency of the vibration. 
The electrical energy loss per cycle Ei oss is 


Eioss — I I(t) 2 (R c + R L )dt 
Jo 

V 2 mp T R c + R l 

2 (R c + Rl) 2 + (L c u>) 2 ’ 


( 7 ) 


where I(t) is the current through the coil and V amp is the amplitude of the voltage generated by the coil. 
To maximize Ei oss , Rl should satisfy 


f T c ca R c (L c uj Rc) 
\ 0 ( L c oj < R c ) 


The maximum electrical energy loss is 


f (TV 2 mp )(AL c uj)- 1 ( L c u>R c ) 

Eioss = \ 

l (TRcV 2 mp ){2^/R 2 c + ( L c w) 2 )-i ( L c u<Rc ) 


( 8 ) 


( 9 ) 


Equations (6) and (8) show that the maximum output power and the maximum loss factor do not 
happen at the same load conditions. To bring the two optimal resistances closer, either R c should be 
reduced or L c needs to be increased. Equation (9) shows that increasing L c will reduce the maximum 
energy loss which is detrimental. On the other hand, the coil resistance R c can be reduced by using 
copper wires with a larger cross section. However, the larger wire size limits the number of turns per 
layer, thus reducing V a m p- This tradeoff will be investigated in future work. 


2.4 Resistive-capacitive load 

Previous studies showed that energy dissipation can be improved by connecting a capacitor in parallel 
with the resistive load (Figure 4(b)). 3 The average power dissipated by the pickup coil is 


Pc = .!. / i 2 COM 

1 Jo 

= \Vg mp \ 2 [( CRlCO ) 2 + 1 }Rc 

2 (R l + R c - LcRlCu 2 ) 2 + ( R L R c C + L c ) 2 lo 2 ' 


( 10 ) 


The average electrical output power is 



2 (Rl + R c 


Rl 

L c RlCuj 2 ) 2 + (R l R c C + L c yu 2 ' 


( 11 ) 
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According to a previous study, 17 the maximum P is 


V 2 

p amp 

= ~8R7’ 


when 


C = 


and Rl = R c 


Rc 


R 2 + (L c u) 2 

The total electrical energy loss per cycle E loss is 
Eloss = (Pc + P)T 

_ \Vam P ?T (R l Cuj) 2 R c + (R c + R l ) 


(R l + R c - L c R l Clo 2 ) 2 + ( RlRcC + L c ) 2 u 2 


( 12 ) 

(13) 


(14) 


which is maximum when C = 1 /(L c uj 2 ) and Rl = +oo. In other words, the maximum energy loss is 
achieved when Rl is disconnected and the L-C circuit has a natural frequency that is the same as the 
driving frequency. The capacitor connected in parallel compensates for the electrical impedance due to 
the coil inductance and improves the loss factor. However, its performance degrades when the driving 
frequency shifts away from the L-C resonance. 


3. NUMERICAL COMPARISON OF GALFENOL AND TERFENOL-D 
3.1 Pickup coil size 

The pickup coil is made of AWG 36 copper wire. The flux density through the rod is assumed to be a 
750 Hz, 0.2 T sinusoidal wave, and the pickup coil is short-circuited. The induced voltage on the pickup 
coil increases with the number of layers until the outer layers become insensitive to the flux density 
change on the rod. On the other hand, the coil resistance increases rapidly with respect to the coil size. 
Hence, an optimal number of layers exists to maximize the average electrical output power of the coil. 
Figure 5 shows that a pickup coil with 29 layers dissipates the maximum amount of energy. 


x 10' 4 



Figure 5. Energy dissipated by the pickup coil per cycle with respect to the number of layers in the pickup coil. 


3.2 Flux path 

The permeability of Terfenol-D is close to that of air, so part of the flux generated by the permanent 
magnets leaks through air. The magnetic held through the Terfenol-D spoke is 


Ht 


4 > 

A r + (no/ ht)Aa 


( 15 ) 
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where $ is the flux generated by the permanent magnet, A r is the cross-sectional area of the magne- 
tostrictive rod, n q is the permeability of air, /, it is the permeability of the Terfenol element, and Aa is 
the cross-sectional area of air. The term ht changes with respect to applied stress, but the variation is 
relatively small (Table 1). Hence, the bias magnetic field through the Terfenol-D rod is nearly constant. 

Galfenol has a relatively high magnetic permeability, so almost all the magnetic flux flows through 
the Galfenol cross section. The magnetic field through the Galfenol rod is 


H g 


<f> 

A r HG ’ 


(16) 


where \ig is the permeability of Galfenol. According to Table 1, [iq has a large variation with respect to 
applied stress, thus the magnetic field Hq changes significantly during operation. Previous research 18,19 
showed that the maximum magnetomechanical coupling is achieved when the magnetic field through the 
magnetostrictive materials is constant. 


To prevent the magnetic field in the Galfenol rod from varying, a flux path is added in parallel to 
the Galfenol rod, as shown in Figure 3(a). The permanent magnets are considered as a constant current 
source $ in the equivalent circuit, as shown in Figure 6. The voltage on the coil V is calculated using 
Faraday’s law: 


v = _ N ^£ = N d ®G dR G dTjt) 
dt dRc dT(t ) dt 1 


(17) 


where N is the number of turns of the pickup coil, as optimized in the previous step, is the flux 
through the Galfenol rod, T(t) is the applied stress, and Rq is the magnetic reluctance of Galfenol. The 
term dRc/dT(t) is determined by Galfenobs material property, while the dT(t)/dt term is determined 
by the input stress wave. The only term in Equation (17) that needs to be optimized is 


dAg _ R c <f> 
dR G ~ ( R c + R g ) 2 ' 


(18) 


When the reluctance of the parallel path, f? c , equals Rq , the output voltage V reaches a maximum. 



Figure 6. Equivalent circuit for the Galfenol system. 

3.3 Purely resistive load 
3.3.1 Terfenol-D 

Figure 7 shows the loss factor r\ and average electrical output power P with respect to the rod’s radius 
R r and the load resistance Rl ■ The maximum output power on the resistive load occurs at around 
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Figure 7. (a) Loss factor with respect to load resistance and rod radius (Terfenol-D); (b) average electrical output 
power with respect to load resistance and rod radius (Terfenol-D). 


Rl = 40 Li and R r = 2.1 mm, and the maximum output power is 0.58 W. When the pickup coil is 
shorted, the device cannot output any energy. For small R ri P is relatively small, since the 280 N 
amplitude input force saturates the material. The inductance of the coil can be estimated by 


_ nN 2 A 

I J /-> 

c L 


(19) 


where /i is the permeability of the core material, L is the length of the coil, and A is the cross-sectional 
area of the coil. When R r is large, the impedance of the coil reduces P. The average electrical output 
power reaches a maximum when Rl = \J + (L c w) 2 . 

The peak 77 occurs at Rl = 55 L2 and R r = 3.9 mm. The maximum 77 is 0.34. Compared to the 
P surface, the maximum 77 occurs at a relatively low load resistance for the same R r . When the loss 
factor reaches its maximum, the average electrical output power is 0.32 W, which is only 53.5% of the 
maximum P. 


3.3.2 Galfenol 

Figure 8 shows the P and 77 surfaces for the Galfenol system. Similar to the Terfenol-D rod, there are 
peaks in both the loss factor surface and the average electrical output power surface. The maximum P 
and 77 are 0.21 W and 0.14, respectively, both of which occurring for a 2 mm return path thickness. The 
peak P and 77 occur for load resistances of Rl = 650 and Rl = 500 O, respectively. Galfenol has a 
high magnetic permeability, which creates a large inductance L c in the equivalent circuit. Hence, large 
Rl is needed to reach the maximum P and 77 , and Galfenol has a relatively low average electrical output 
power and loss factor compared to Terfenol-D. 

3.4 Resistive-capacitive load 

This study considers the case when Rl = +00 (Figure 4(b)) and Cl is varied. To reduce the compu- 
tational effort, only a Terfenol-D rod with a 3.9 mm radius and a Galfenol rod with a 2.5 mm radius 
are modeled. Figure 9 shows that the loss factor can be improved by impedance matching. The maxi- 
mum loss factors for Galfenol and Terfenol-D at resonance are 0.59 and 1.83, respectively. But this L-C 
resonant electric circuit tends to create a narrow loss factor peak in the frequency domain. 
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Figure 8. (a) Loss factor with respect to load resistance and rod radius (Galfenol); (b) average electrical output 
power with respect to load resistance and rod radius (Galfenol). 
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Figure 9. (a) Loss factor of the Galfenol rod (2.5 mm radius) with respect to varying capacitance; (b) loss factor 
of the Terfenol-D rod (3.9 mm radius) with respect to varying capacitance. 
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4. SUMMARY AND FUTURE WORK 


This study provided a modeling framework and design insight for developing a magnetostrictive vibra- 
tion ring. This device would be mounted within a driveline and transfer mechanical vibration energy 
to electrical energy to create a damping effect. The work presented focused on optimizing the mag- 
netostrictive elements, which are either Terfenol-D or Galfenol. A parametric study considering loss 
factor and average electrical output power was conducted using a 2D axisymmetric finite element model. 
A 750 Hz, 280 N amplitude excitation force was applied to both rods. For purely resistive loads, the 
maximum loss factors of Galfenol and Terfenol-D are 0.14 and 0.34, respectively. The maximum average 
electrical output power P for Galfenol and Terfenol-D is 0.21 W and 0.58 W, respectively. The difference 
is attributed to the relatively high permeability of Galfenol. This leads to an increased inductance which 
restricts current variation in the electric circuit. By using an L-C resonant circuit, the maximum loss 
factors for Galfenol and Terfenol-D increase to 0.59 and 1.83, respectively. 

In the next stage of this work, an experiment will be conducted to validate the modeling results. 
While this work focused on a single frequency, future modeling efforts will investigate the bandwidth of 
the system. Results from the present study predict that P and rj could be simultaneously maximized by 
using a copper coil with larger cross section. This improvement will be incorporated in future designs. 

ACKNOWLEDGMENTS 

We wish to acknowledge the member organizations of the Smart Vehicle Concepts Center, a National 
Science Foundation Industry/University Cooperative Research Center (www.SmartVehicleCenter.org) 
established under NSF Grant IIP-1238286. 

REFERENCES 

1. V. M. Asnani, T. L. Krantz, D. C. Delap, and D. B. Stringer, “The Vibration Ring Seedling Fund 
Phase 1 Final Report,” NASA TM-218337, 2014. 

2. D. Davino, A. Giustiniani, C. Visone, and A. Adly, “Experimental analysis of vibrations damping 
due to magnetostrictive based energy harvesting,” Journal of Applied Physics 109(7), p. 07E509, 
2011. 

3. J. Yoo, A. Murray, and A. B. Flatau, “Evaluation of magnetostrictive shunt damper performance 
using iron (Fe)-gallium (Ga) alloy,” in SPIE Smart Structures and Materials & Nondestructive Eval- 
uation and Health Monitoring , pp. 90573I-90573I, International Society for Optics and Photonics, 
2014. 

4. T. Ueno and S. Yamada, “Performance of energy harvester using iron-gallium alloy in free vibra- 
tion,” Magnetics, IEEE Transactions on 47(10), pp. 2407-2409, 2011. 

5. J. H. Yoo and A. B. Flatau, “A bending-mode Galfenol electric power harvester,” Journal of Intel- 
ligent Material Systems and Structures 23(6), pp. 647-654, 2012. 

6. V. Berbyuk, “Vibration energy harvesting using Galfenol-based transducer,” in SPIE Smart Struc- 
tures and Materials & Nondestructive Evaluation and Health Monitoring , pp. 86881F-86881F, In- 
ternational Society for Optics and Photonics, 2013. 

7. L. Wang and F. Yuan, “Vibration energy harvesting by magnetostrictive material,” Smart Materials 
and Structures 17(4), p. 045009, 2008. 

8. M. B. Moffett, A. E. Clark, M. Wun-Fogle, J. Linberg, J. P. Teter, and E. A. McLaughlin, “Char- 
acterization of Terfenol-D for magnetostrictive transducers,” Journal of the Acoustical Society of 
America 89(3), pp. 1448-1455, 1991. 

9. A. E. Clark, M. Wun-Fogle, J. B. Restorff, and T. A. Lograsso, “Magnetostrictive properties of 
Galfenol alloys under compressive stress,” Materials Transactions 43(5), pp. 881-886, 2002. 


Proc. of SPIE Vol. 9433 94330C-10 


Downloaded From: http://proceedings.spiedigitallibrary.org/ on 04/07/2015 Terms of Use: http://spiedl.org/terms 


10. R. Kellogg, A. B. Flatau, A. Clark, M. Wun-Fogle, and T. A. Lograsso, “Temperature and stress 
dependencies of the magnetic and magnetostrictive properties of Feo. 8 iGao. 19 ,” Journal of Applied 
Physics 91(10), pp. 7821-7823, 2002. 

11. M. Wun-Fogle, J. Restorff, and A. Clark, “Magnetomechanical coupling in stress-annealed Fe-Ga 
(Galfenol) alloys,” Magnetics, IEEE Transactions on 42(10), pp. 3120-3122, 2006. 

12. http:/ /www. etrema-usa.com/documents/terfenol. pdf 

13. S. Chakrabarti, Modeling of 3D magnetostrictive systems with application to Galfenol and Terfenol- 
D transducers. PhD thesis, The Ohio State University, 2011. 

14. Z. Deng and M. J. Dapino, “Characterization and finite element modeling of galfenol minor flux 
density loops,” Journal of Intelligent Material Systems and Structures , p. 1045389X14521703, 2014. 

15. E. E. Ungar and E. M. Kerwin Jr, “Loss factors of viscoelastic systems in terms of energy concepts,” 
Journal of the Acoustical Society of America 34(7), pp. 954-957, 1962. 

16. ASTM D5992-96(2011), “Standard guide for dynamic testing of vulcanized rubber and rubber-like 
materials using vibratory methods,” ASTM International, West Conshohocken, PA, 2011. 

17. Z. Deng and M. J. Dapino, “Modeling and design of Galfenol unimorph energy harvester,” in SPIE 
Smart Structures and Materials & Nondestructive Evaluation and Health Monitoring , pp. 90572A 
90572A, International Society for Optics and Photonics, 2014. 

18. L. Weng, T. Walker, Z. Deng, M. Dapino, and B. Wang, “Major and minor stress-magnetization 
loops in textured polycrystalline Fegi. 6 Ga 18 4 Galfenol,” Journal of Applied Physics 113(2), 
p. 024508, 2013. 

19. A. Mahadevan, P. Evans, and M. Dapino, “Dependence of magnetic susceptibility on stress in tex- 
tured polycrystalline Fe 8 i. 6 Gai 8.4 and Fe 7 g.iGa 2 i .9 Galfenol alloys,” Applied Physics Letters 96(1), 
pp. 012502-012502, 2010. 


Proc. of SPIE Vol. 9433 94330C-11 


Downloaded From: http://proceedings.spiedigitallibrary.org/ on 04/07/2015 Terms of Use: http://spiedl.org/terms 


